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Abstract. We discuss the first passage time problem in the semi-infinite 
interval, for homogeneous stochastic Markov processes with Levy stable jump 
length distributions A(x) ~ ^ £), namely. Levy flights (LPs). In 

particular, we demonstrate that the method of images leads to a result, which 
violates a theorem due to Sparre Andersen, according to which an arbitrary 
continuous and symmetric jump length distribution produces a first passage time 
density (FPTD) governed by the universal long-time decay ~ t~^^^. Conversely, 
we show that for LPs the direct definition known from Gaussian processes in 
fact defines the probability density of first arrival, which for LPs differs from the 
PPTD. Our findings are corroborated by numerical results. 
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Levy flights (LFs) and Levy walks (LWs) are the prime example in the 
investigation of non-standard transport processes whose stationary solution docs not 
converge towards the Boltzmann form ^EIEIISI- Being subject to the generalised 
central limit theorem E| , LFs correspond to a Markov process in which extremely 
long excursions can occur with appreciable probability, whereas in LWs long excursions 
are penalised through a time cost introduced via a spatiotemporal coupling ^7, jS]. 
Applications of LFs and LWs range from the famed flight of an albatross the 
spreading of spider- monkeys (TD| . or the grazing patterns of bacteria over 
economical data 12 to molecular collisions ^31 and plasmas ^3]. Despite their 
broad usage, the detailed behaviour of even the simpler, uncoupled LF processes, 
on which we concentrate in the following, in external potentials and under non- 
trivial boundary conditions is still not fully explored. Thus, there have recently been 
discovered bifurcations between multimodal states of the probability density function 
(PDF) of LFs in steeper than harmonic external fields, in whose presence also the 
variance becomes finite [151 , and rich band structures have been reported for LFs 
in periodic potentials [TTj . 

Of particular interest in random processes is the first passage time density (FPTD) 
^1^11201 IS]- For LFs, the FPTD was determined through the method of images 
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on a finite domain in reference 22 , and with similar methods in reference j23| . These 
methods lead to results for the FPTD in the semi-infinite domain, whose long-time 
behaviour explicitly depends on the Levy index a. In contrast, a theorem due to 
Sparre Andersen proves that for any discrete-time random walk process starting at 
Xq ^ with each step chosen from a continuous, symmetric but otherwise arbitrary 
distribution, the FPTD asymptotically decays as ~ n~^^^ with the number n of steps 
1^1^1231, being fully independent of the index of the LP, i.e., universal. In the case 
of a Markov process, the continuous time analogue of the Sparre Andersen result reads 

pit)-t-'^\ (1) 

The analogous universality was proved by Frisch and Prisch for the special case in 
which an absorbing boundary is placed at the location of the starting point of the 
LP at f > 1201 1 and numerically corroborated by Klafter and Zumofen |27]. In the 
following, we demonstrate that the method of images is generally inconsistent with the 
universality of the PPTD, and therefore cannot be applied to solve PPTD-problems 
for LPs. We also show that for LPs the PPTD differs from the PDP for first arrival. 

Let us start by recalling that an unbiased LP can be defined through the space- 
fractional diffusion equation for the PDP W{x,t) pi l^l^ 

where we define the fractional derivative d" /d\x\°' by its Fourier transform. (Here and 
in the following, we restrict ourselves to 1 < a < 2.) In position space, the fractional 
derivative is defined in terms of the convolution (see jl5j for the case a = 1) 

d" , D [°° W{x',t) 



7ra 
cos — 
2 



n dx^ J_^\x-x'\°'-^ --'^ " ™ 2 '^^^ 

Equivalently, LPs can be described in terms of continuous time random walks with 
long-tailed jump length distributions A(a;) ~ /\x\^^°' [31301. The associated PDF 
W{x,t) for natural boundary conditions (lim|2,|^oo W{x^t) = 0) with initial condition 
5{x) is the Levy stable law W{x,t) = exp (-ifcx - i:>|fc|"i) dfc/(27r) (Sj 0. In 

Fourier-Laplace space [S], this PDF corresponds to W{k,s) = (s-|-D|fc|") ^. A 
characteristic of LPs is the divergence of the variance of both W{x,t) and \{x). 
Equipping equation (j^J with a i5-sink of strength pfa(i), we obtain the diffusion-reaction 
equation for the non-normalised density function f(x,t), 

^/(x, t) = D^fix, t) - pUm^l (4) 

from which by integration over the mirestricted space, we find the quantity 

Pfa(i) = -^ j /(x,i)dx, (5) 

i.e., Pfa(i) is the negative time derivative of the survival probability. By definition 
of the sink term, pfa(i) is the PDF of first arrival: once a random walker arrives 
at the sink, it is annihilated. By solving equation through standard methods 
(determining the homogeneous and inhomogeneous solutions) , it is straightforward to 
calculate the solution / in terms of the propagator W of equation Q with initial 
condition f{x,Q) — 5[x — xq) yielding f{k,u) — [e'^^'^o + p{a{u)] / {s + D\k\"), from 
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which, in turn, we find that pia{t) satisfies the chain rule (pfa impficitly depending on 

W{-xo,t) = f pUt)W{0, t - T)dr (6) 
Jo 

which corresponds to the Laplace space relation pfa(s) = W{—xq,s)/W{Q,s). 
Equation (O is well-known and for any sufficiently well-behaved continuum diffusion 
process is commonly employed to define the FPTD OT] . 

For Gaussian processes with propagator W{x,t) = l/v47ri5texp (— a:^/[4£>t]), 
one obtains by direct integration of the difi^usion equation with appropriate boundary 
condition the FPTD [H] 

p{t) = xoiAnDt^)-'/^ exp {-xl/{ADt)) , (7) 

including the asymptotic behaviour p{t) ^ t~^^^ for t ^ Xq/{4D). In this Gaussian 
case, the quantity pf^it) is equivalent to the FPTD. From a random walk perspective, 
this is due to the fact that individual steps are of the same increment, and the jump 
length statistics therefore ensures that the walker cannot hop across the sink in a 
long jump without actually hitting the sink and being absorbed. This behaviour 
becomes drastically different for Levy jump length statistics: There, the particle 
can easily cross the sink in a long jump. Thus, before eventually being absorbed, 
it can pass by the sink location numerous times, and therefore the statistics of 
the first arrival will be different from the one of the first passage. In fact, with 
W{x, s) = (27r)-i e"^^ (s + i:>|fc|")"^ dfc, we find 

/^(l-COsfcxo)/(5 + Dfc")dfc 

P'^^'^ = ' /-l/(.+m^)dfc — 

by use of the de Moivre identity exp(iz) = cosz -I- isinz. With /q°°(s + Dk")^^dk = 
ws^/"-^/{aD^/°'sm{TT/a)) and 

1 - cos kxo r((2 - a) sin {tt(2 - a) 12) xf^^ 

^ ; 7-;^ — ■ s^O,a>l, 

s + Dk°' {a-l)D 

we obtain the limiting form 

Pfa(s) ~ 1 - x^isi-i/"i?-i+i/"A(a), (9) 

where A(a) = aT{2 — Q;)sin(7r(2 — q;)/2) sin(7r/Q!)/(Q; — 1). We note that the same 
result is obtained by the exact expressions for W{xo,s) and 14^(0, s) in terms of Fox 
H-functions and systematic expansion The inverse Laplace transform of the 

small s-behaviour @ can be obtained by completing ^ to an exponential, and then 
computing the Laplace inversion by the identification — Hq'^[z\{0, 1)] with the Fox 
i?-function , for which the exact Laplace inversion can be performed [33] . Finally, 
a series expansion of this result leads to the long-t form 

PfaW-C(a) ^,."Vl/. ^ (10) 

with C{a) = ar(2-a)r(2- 1/a) sin(7r[2 - a]/2) sin2(7r/a)/(7r2(a - 1)). Clearly, in 
the Gaussian limit, the required asymptotic form p{t) ^ xo/v4^7TDt^ for the FPTD 
is consistently recovered, whereas in the general case the result H1Q(I is slower than 
in the universal FPTD behaviour in equation (^), as it should as the (5-trap used in 
equation Q to define the first arrival for LFs is weaker than the absorbing wall used 



oo 



Letter to the Editor 



4 



a=1.2, x„=0.0, w=0.3 
a=1.2, Xo=0.3, w=0.3 
a=1.2, x„=1.0, w=0.3 
a=1.2, Xq=10.0, w=1.0 
a=1.8, Xn=10.0, w=0.25 
tt! ° a=1.2 



□ □ 



igt 



Figure 1. First arrival PDF for a = 1.2 demonstrating the (-2+i/a! scaling, 
for optimal trap width w = 0.3. For comparison, we show the same scaling for 
a = 1.8, and the power-law t~^^'^ corresponding to the FPTD. The behaviour for 
too large w = 1.0 shows a shift of the decay towards the —3/2 slope. Note that 
on the abscissa we plot Igtp(t). Note also that for the initial condition xq = 0.0, 
the trap becomes activated after the first step, consistent with 127) . 



to properly define the FPTD. For LFs, the PDF for first arrival thus scales like (jlO|l 
(i.e., it explicitly depends on the index a of the underlying Levy process), and, as 
shown below, it differs from the corresponding FPTD. 

Before calculating this FPTD, we first demonstrate the validity of equation (|10|l 
by means of a simulation the results of which are shown in figure ^ Random jumps 
with LF jump length statistics are performed, and a particle is removed when it hits a 
certain interval of width w around the sink; for our simulations we found an optimum 
value w ~ 0.3. As seen in figure H (note that we plot \gtp{t)\) and for analogous 
results not shown here, relation (|10|l is nicely fulfilled for all 1 < a < 2, whereas for 
larger w, the slope increases. 

The proper dynamical formulation of an LF on the semi-infinite interval with an 
absorbing boundary condition at a; = 0, and thereby the determination of the FPTD, 
has to make sure that in terms of above random walk picture jumps across the sink 
are forbidden. This can be consistently achieved by setting f{x,t) = on the left 
semi-axis, i.e., actually removing the particle when it crosses the point a; = 0. This 
formally corresponds to the modified dynamical equation 



dt K dx^ Jq Ix-x'I"-! dx^' 
in which the fractional integral is truncated to the semi-infinite interval. After Laplace 
transformation and integrating over x twice, one obtains 

K{x — x', s)f{x', s)dx' — {x — xi:,)Q{x — X()) — xp{s) — T{0, s), (12) 







where p{t) is the FPTD and the kernel K{x^s) = sxQ{x) — {k\x\'^ ^). This equation 
is formally of the Wiener- Hopf type of the first kind |34| . After some manipulations 
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x_0=0.10 + 
x_0=1 .00 X 
x_0=10.0 X 
x_0=100.0 □ 

t"(-1.5) 

f(-1-5) - — 
t**(-1-5) - 
t"(-1-5) 
First arrival 
■^■□ Images method 




Figure 2. Numerical results for the FPTD process on the semi-infinite domain, 
for an LF with Levy index a = 1.2. Note that on the abscissa, we plot tp{t). For 
all initial conditions xq = 0.10 1.00, 10.0, and 100.0 the universal slope —3/2 in 
the logj^Q-logjQ plot is nicely reproduced, and it is significantly apart from the two 
slopes predicted by the images method and the direct definition of the FPTD. 



similar to those applied in reference |23, we arrive at the asymptotic expression 
p(s) ~ 1 — Cs^/^, where C = const, in accordance with the universal behaviour 
and with the findings in reference Thus, the dynamic equation 1)11(1 consistently 
phrases the FPTD problem for LFs. We note that due to the truncation of the 
fractional integral it was not possible to modify the well-established Griinwald- 
Letnikov scheme [SS] to numerically solve equation (|llll with enough computational 
efficiency to numercially obtain the direct solution for f{x, t). However, to corroborate 
the validity of the Sparre Anderson-universality, we perform a simulation of an LF in 
the presence of an absorbing wall, i.e., random jumps with LF jump length statistics 
are performed along the right semi-axis, and a particle is removed when it jumps 
across the origin to the left semi-axis. Results of such a detailed random walk study 
are displayed in figures|lland|31 The expected universal t~^/^ scaling is nicely obtained 
for various initial conditions and a. Clearly, the scaling for the first arrival as well as 
the image method-FPTD derived below are significantly different. 

We now demonstrate that the method of images produces a result, which is neither 
consistent with the universal behaviour of the FPTD nor with the behaviour of the 
PDF of first arrival H10() . Given the initial condition S{x — xq), the solution fini{x,t) 
for the absorbing boundary value problem with the analogous Dirichlet condition 
/im(0,t) — according to the method of images is given via the difference |2()l I21| 

fi^ix, t) = W{x - xo, t) - W{x + xo, t), (13) 

in terms of the free propagator W , i.e., a negative image solution originating at 
—Xq balances the probability flux across the absorbing boundary. The corresponding 
pseudo-FPTD is then calculated in the same way as in equation Q. For the image 
solution in Fourier-Laplace space, we obtain 

fUk,s) = [2isin(fca;o)]/(s + ^|/cr), (14) 
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Figure 3. Same as in figure 1^ for a = 2.0, 1.5, 1.0, and 0.6, and for the initial 
condition xq = 10.0. Again, the universal ~ t^^/^ behaviour is obtained. 



for a process which starts at > and takes place in the right half space. In Laplace 
space, the image method-FPTD becomes 



Pim(s) = 1 - 



dx 



dk 
— f 
27r 



2i sin kxQ 
s + D\k\°'' 



(15) 



After some transformations, we arrive at 

Pi,„(s) = 1 - 2/7r /;° sin {^s^/'^xo/ D^/-) / (1 + C)] ■ (16) 

In the limit of small s, this expression reduces to pim(s) ~ 1 — A{a)xoD~^^°'s^^°' , with 
A(a) = (2/7r) {I + d^ = 2/(asin(7r/a)). Following the same procedure as 

outlined above, we find the long-t form 

Pi^{t) - 2r(l/a)a;o/ 



(17) 

2, expression (fT7|l produces 



naD^/''t^+^/°' 

for the image method-FPTD. In the Gaussian limit a 
Pim(i) ~ xo/v47rDt^, in accordance to equation Conversely, for general 1 < a < 2, 
Pim(i) according to equation (fT7|l would decay faster than ~ t~^^'^. The method of 
images breaks down for LFs due to their special non-local nature, displayed by the 
integrals in equations (0) and lO, and Hll() . namely having a long-tailed jump length 
distribution. This leads to leapovers beyond the absorbing boundary. The method of 
images is expected to work when the boundary is also a turning poing of the trajectory, 
as actually happens for nearest neighbours random walks, or the Wiener process. 

Qualitatively, the following argument may be brought forth in favour of the 
observed universality of the LF-FPTD: the long-time decay is expected to be governed 
by short-distance jump events, corresponding to the central region of very small jump 
lengths for the Levy stable jump length distribution. But this region is, apart from 
a prefactor, indistinguishable from the Gaussian distribution, and therefore the long- 
time behaviour should in fact be the same for any continuous jump length distribution 
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X{x). In fact, the universal law JQ) can only be modified in the presence of non- 
Markov effects such as broad waiting time processes or spatiotemporally coupled walks 
121 El U\ Uni inn ■ In terms of the special case covered by the theorem of Frisch and 
Frisch |26|. in which the absorbing boundary coincides with the initial condition, we 
can understand the general situation for finite xq > 0, as in the long-time limit, the 
distance Xo becomes negligible in comparison to the diffusion length (|a;(t)|) ~ t^^": 
therefore the asymptotic behaviour is necessarily governed by the same universality. 

Concluding, we demonstrated that the method of images, which has been 
developed a powerful tool in Gaussian diffusion also beyond the homogeneous case 
|20l 1^ and in the presence of long-tailed waiting times 12 [211 EZj, fails for LF 
processes, leading to a false result for the FPTD. Moreover, we showed that for such 
broad jump length statistics, the PDF of first arrival at a point differs from the FPTD. 
We also provided a framework in terms of a truncated fractional diffusion equation to 
solve the FPTD problem for an LF. This study is expected to significantly contribute 
to the understanding of the, at instances, non-trivial behaviour of LFs. 
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